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For a Bose condensate in a double- well potential or with two Josephson-coupled internal states, the 
condensate wavefunction is a superposition. Here we consider coupling two such Bose condensates, 
and suggest the existence of a joint condensate wavefunction, which is in general a superposition 
of all products of the bases condensate wavefunctions of the two condensates. The corresponding 
many-body state is a product of such superposed wavefunctions, with appropriate symmetrization. 
These states may be potentially useful for quantum computation. There may be robustness and 
stability due to macroscopic occupation of a same single particle state. The nonlinearity of the 
.. ■ condensate wavefunction due to particle-particle interaction may be utilized to realize nonlinear 

^.«: I quantum computation, which was suggested to be capable of solving NP-complete problems. 

^ . PACS numbers: 03.75.Fi, 05.30.-d, 03.67.Lx, 03.65.Ud 

The Hilbert space of the quantum state of a composite system of coupled particles is the tensor product of the 

QQ Hilbert spaces of these particles. That is to say, a quantum state of a composite system is in general a superposition 

K*" . of all the products of the bases states of the constituent particles. This leads to the situation that the state of 

^^ ' a composite system may not be factorized to be a direct product of the states of the constituent particles. This 

L.^ , so-called quantum entanglement has been a central attention in foundations of quantum mechanics ever since the 

^— V ' early days pKB- Moreover, in recent years, it was found that superposition and entanglement of quantum states give 

^-H . rise to powerful quantum information processing. For example, a quantum computer is much more efficient than a 

0^ ' classical computer in solving certain problems such as factoring large numbers |B|J6|. It is, however, highly challenging 

^^ ] to physically build a quantum computer, a major obstacle being decoherence, i.e. the fragility of a superposed or 

C^ ' entangled state towards coupling with the environment, although the quantum error-correction code diminishes this 

Oh! difficulty [[7|. Stimulated by both its power and the challenge in its physical implementations, the field of quantum 

■4.^ ' information is attracting more and more attention. 

Q Here we extend the consideration to coupled many-particle systems, each of which is composed of identical particles, 

r^ ' but the coupled systems are distinguishable with each other. Specifically, we consider coupled Bose condensates, each 
qh] of which is in a double- well potential or has two Josephson-coupled internal states (note that we refer to Josephson- 
IL* '■ coupled condensates as one condensate, since in this case, all the particles, even though from different sources originally, 
. ^ \ are overlapping identical particles Q). Because of coupling between each particle in one condensate and each particle in 
j^ ' the other condensate, we consider the mean-field many-body state which is a product of copies of a same wavefunction 
Jh \ of inter-condensate pair of particles, with symmetrization. With the flexibility of manipulation, this situation may be 
- - ' realizable in the context of Bose-Einstein condensation of trapped atoms ||9|-[l^ . 

A Bose condensate is well described by the condensate wavefunction, which, in a mean field theory, is the single 
particle wavefunction in which the condensation occurs. It can also be defined as the spontaneous gauge symmetry 
breaking (SGSB) of the field operator, or through the off-diagonal long-range order (ODLRO) of the one-particle 
reduced density matrix. It is a superposition for a condensate in a double-well potential or with two Josephson- 
coupled internal states. In our consideration of the coupled condensates, the wavefunction of each inter-condensate 
pair of particles becomes the joint condensate wavefunction, which can entangle the condensate wavefunctions of 
individual condensates. 

With superposition and entanglement, it may be possible to use condensate wavefunction to implement a qubit in 
quantum computation. The first possible advantage is the intrinsic robustness and stability. It may be viewed as a 
natural realization of a fault tolerance prescription based on symmetrization p3| . To have the best accuracy, each 
Bose system had better be non-interacting. In principle, this is possible since the interaction can be tuned in the 
case of atomic condensate El. On the other hand, it was noted that incorporating nonlinear evolution in quantum 
computation leads to computational power of solving NP-complete and #P problems in polynomial time ||l5| . This 
result has remained as an academic curiosity, since quantum mechanics is linear. Now, for a weakly interacting 
Bose system, the condensate wavefunction satisfies a nonlinear Schrodinger equation, often called Gross-Pitaevskii 



equation ||I^. Thus Bose condensates may make it possible to implement nonlinear quantum computation, for the 
sake of special power. Therefore Bose statistics may be a resource of both fault tolerance and computational power. 
We will illustrate the ideas in terms of Bose-Einstein condensation of trapped atoms, mainly using a condensate in a 
double-well potential, with the two bits represented by the condensate wavefunctions localized at the two wells. We 
also discuss a condensate with two Josephson-coupled internal states, which might encode the bits. But there are 
difficulties due to coupling with the motional degree of freedom and the nonlinearity. However, the physical problems 
are still interesting. 

The paper is organized in the following way. Section 2 is an overview of the concept of condensate wavefunction. 
In Appendix A, we give a derivation of Gross-Pitaevskii equation from ODLRO. Bose-Einstein condensation in a 
double-well potential is a prototype in discussions from Sees. 3 to 5. As a prelude, superposition of condensate 
wavefunction is discussed in Subsec. 3.1. In Subsec. 3.2, we consider coupling two different condensates, and discuss 
possible entanglement between them. Many-body Hamiltonian and the equation of motion of the joint condensate 
wavefunction are given in Sec. 4. Section 5 contains illustrative schemes of one-bit and two-bit operations in terms of 
condensates in double-well potentials. In Sec. 6, we discuss, in parallel to Sees. 3 to 5, the case of spinor condensate, 
which has two components with different internal states. We suggest that coupling two spinor condensates gives rise 
to a four-component condensate wavefunction. Section 7 contains some additional remarks. A summary is given in 
Sec. 8. 

2. INTRODUCING CONDENSATE WAVEFUNCTION 

In the following, we review three approaches to the condensate wavefunction. 

(i) In the mean field theory, which becomes exact in the absence of the interaction, the Bose-condensed state is in 
a Hatree form, 

^(ri,---,riv) = (/'(ri)---<?!)(riv), (1) 

where (j){ri) is the single particle state, which turns out to be the condensate wavefunction. One can define 

$(r) = VN(t){r). (2) 

(ii) In the approach of SGSB (l7|^, 

$(r) = (^(r)), (3) 

where ip{r) is the boson field operator. The particle number is not conserved. One may use (0) to define (j>{r), with 
N now being mean particle number. 

(iii) The general criterion for Bose-Einstein condensation is the ODLRO of the one-particle reduced density ma- 
trix @,|l|, 

(r'Ipilr) = (V3(r')^t(^)) ^ Tr[pi^{r')^^ (r)], (4) 

where p is the density matrix. ODLRO can be expressed as 

(r'|pi|r) ^$(r')$(r)* 7^0 (5) 

as \r — r'\ -^ oo. This applies to both particle number conservation and non-conservation cases. In case of particle 
number conservation, $(r) = \'^^o{r), where <I'o('") is the eigenfunction of pi with the largest eigenvalue A, of the 
order o{N). In case of particle number non-conservation, ^(r) = (-(/'(r)). When N — + cx3, there is no practical 
difference in using the particle number non-conserved state or the particle number conserved state. 
The many-body Hamiltonian is 

N 



iI = ^Mr,)+^C/(r,-r,), (6) 

where 



i<j 



kr)^-^ + V{r) (7) 

is the single particle Hamiltonian, U{ri — Tj) is the particle-particle interaction, to is the particle mass, V{r) is the 
external potential, e.g. the trapping potential in case of the trapped atoms. In terms of the field operator, the 
many-body Hamiltonian is 

H = / drijj\r)h{r)i;{r) + ^ dridr2V'^(ri)V'^(r2)t/(ri - r2)^(r2)^(ri), (8) 

which leads to the equation of motion of ^/'(r), 

z/i^4^ = [--^^ + V^(r) + f^j\r',t)U{r--r')ti;{r',t)dr'U{r,t). (9) 

ot 2m J 

With weak interaction, one can use s-wave approximation U{r — r') — gS{r — r'), with g — 4Trh?ri/m, where i] is 
the s-wave scattering length. Under the SGSB ansatz, one replaces ^p{r,t) with y/N(f> to obtain the Gross-Pitaevskii 
Equation for the condensate wavefunction. 



It was shown that the ground state energy and density given by Gross-Pitaevskii equation become exact as the particle 
number N tends to be infinity while Nrj is fixed p2], and the error is about 1% under current experimental condition 

Gross-Pitaevskii equation can also be obtained by using the Hatree wavefunction {M in the many-body Schrodinger 
equation, with the Hamiltonian (o) |16| . It can also be obtained phenomenologically from a Ginzburg-Landau func- 
tional uW. In Appendix A, we give a derivation from ODLRO. 

3. SUPERPOSITION AND ENTANGLEMENT OF CONDENSATE WAVEFUNCTIONS 
3.1. Superposition of condensate wavefunctions 

In discussing Josephson-like effect, e.g. a condensate in a symmetric double-well potential, a widely used ansatz 
is that the total condensate wavefunction is a superposition of two bases wavefunctions, with time-dependent coeffi- 
cients |gj,g5[. As a prelude to the next subsection, here we make some discussions based on the construction of the 
field operator and then imposing SGSB or ODLRO. 

The basis set of the single particle wavefunctions is {4>a,n}, where a denotes the energy level, n = 0, 1 denotes 
the two wells. Because of finiteness of the barrier, | J 0* g(r)(/>Q_i(r)(irp = e << 1, i.e. the two states are nearly 
orthogonal though not exactly so. The field operator can be constructed as 

n— 0,1 Q 

where aa,n{t) is the annihilation operator corresponding to the single particle state (pa.ni'f')- Therefore 

i>{r,t)^Mr,t)+Mr,t), (12) 

where 

'4'n{r,t) ^^(j)a,„ir)daAt), (13) 



with n = 0, 1. By making SGSB average of (p_2[), we have 

$(r,t) = $o(r,i) + *i(r-,0, (14) 

with ^{r,t) = {ip(r,t)), ^nif^t) = {^pnif^t)). Therefore, a general condensate wavefunction is a superposition of the 
condensate wavefunctions corresponding to the two wells |2q] . 



A justification can also be made in terms of ODLRO. With (12), tfie one-particle reduced density matrix is 

(r'Ipilr) = (r', 0|pi|r, 0) + {r', 0|pi|r, 1) + (r', l|pi|r, 0) + (r', l|pi|r, 1), (15) 

where {r',n'\pi\r,n) = {ipn'{r' ,t)^l^{r,t)), with n ^ 0,1, n' ~ 0, 1. The existence of ODLRO implies 

(r',n'|pi|r,n)^<&„,(r',t)$:(r,t). 

Hence 

(r'|pi|T-)^(<5'o(r-',t) + <5'i(r',i))($o(r,<) + $i(r,t))*, 

which indicates the existence of the superposed condensate wavefunction $(r, t) = ^o{r, t) + <i>i(r, t). 

In a mean field theory, each particle occupies the single particle ground state, hence the condensate wavefunction 
is given by the superposition of the two single particle ground states at the two wells. This is called two-mode 
approximation in p4| . So 

,^(r,i) = ^^=co(t)0ao,o(r-)+ci(i)<^ao,iW, (16) 



where (/)Q,o,„(r) is the single particle ground state at well n. The many-body ground state is M| 

*(ri, • • • , rjv, i) = </>(ri, i) • • • (l){rN, t). (17) 

For brevity, we write (|l6| ) as 

0(r,t) = co(i)<^o(r)+ci(i)0i(r). (18) 

3.2. Entanglement of condensate wavefunctions 

Consider two coupled Bose systems a and 6, say, in a double- well potential (Fig. 1). For simplicity but without 
lose of generality, suppose each consists of N identical particles. For system a, the wells are denoted as n", while for 
system 6, the wells are denoted as n''. These two Bose systems should be non-overlapping if they are composed of a 
same kind of particles, but can be mixed if they are composed of two different kinds of particles. 

The existence of a joint order parameter can be seen in general by considering the product of the two field operators. 

^(r^ t)4^{r\ t)^J2J2 ^"" (^"' *)^"' (^^ ^)' (19) 

Imposing SGSB, one has 

$(r^ r^ t) = E E *«",«' (^"' ^'' ^)' (20) 

where ^{r°-,r'',t) = (^p{r°-,t)iil;{r^,t)), $„a „b(r°,r'',<) = {iljn'^{r°-,t)iil;„b{r'',t)). Because of coupling between a and 
6, in general the many-body state cannot be factorized to be a direct product of that of a and that of b, and there is 
interaction energy between a and b. Consequently, ^{r°',r^,t) cannot be factorized to be a product of a condensate 
wavefunction of system a and a condensate wavefunction of system b. 

The justification can also be made in terms of ODLRO. One can define a one-particle-pair reduced density matrix, 

= E E E E(^"'' ^^"'' "-"'^ "''|pl|r^ n^ r^ n"), (21) 

where {r'" ,n'" ,r'>' ,n'>'\pi\r'',n'',r'',n'') = (^n-'ir"' ,t)i;^,>{r''' ,t)^l,{r\t)^la{r'',t)). With ODLRO, 

ii^n^-ir'^', i)Vi„- (r'', t)ipl (r^ t)i^l (r^ t)) ^ $„„,,„., (r-^', r"' , i)<„ „, {r\r\ t). (22) 



Consequently, (^ approaches 

implying that there exists a joint condensate wavefunction J2n'^ n^ '^n'^,n>'{i''' ii'^ ^i)- 
For convenience, we may write (pO[) as 



0(r^r^t) = EE'^«^"''(^"'^^^)' (23) 

with </.(r^r^i) = *(r^r^i)/7V, 0„.,„.(r", r^ t) = $„„,„. (r", r^ <)/iV. 
Consider the following Bose-condensed state, 

*K, ..., r^; rj, ...,<; t) = -^ E['^"S,"!;(^i ' ^^i' *) ' ' ' '^"S.^S (^^, ^W, *)], (24) 

V JV. p 

where 

0„g,„.(r^r^t) = EE^«"«^W'^"8-"(^")'^og,„^(r''). (25) 

The symmetrization in state (E4h represents all possible ways of pairing. The systems Bose-condense into a coupled 
single particle pair state, and (|25|) can be used as the joint condensate wavefunction. In this state, particles in a 
same system occupy a same ground state, while maintain coupling with particles in the other system. It is a natural 
generalization of the ansatz dT^). However, it remains to be rigorously examined whether this state is indeed the 
many-body ground state. 

Like using two-mode approximation for a single condensate, one may make two-mode approximation for each 
condensate, and consider the coupling between them. Then, in consistent with the above mean-field inany-body 
state, the joint condensate wavefunction is 

For brevity, we write 

</.(r^r^i) = EE^""»^(*)'^""(^")'^»^(^'')■ (27) 



4. MANY-BODY HAMILTONIAN AND EQUATION OF MOTION 

For two coupled Bose systems a and 6, the general form of the Hamiltonian is 

N N 

H = Y. h{r-) + E Ui< - ^") + E ^(^^) + E ^(^fe - ^') + E ^(^" - ^fe)' (28) 

i—1 i<j k—1 k<l i,k 

where /i(rf ) and h{r\) are single particle Hamiltonians of a particle in a and h, respectively, as given in (0). U{rf — r'^) 
is the particle-particle interaction within a, while U{r^^ — r^) is the interaction within b. W{rf — r^,) is the interaction 
between a particle in a and a particle in b. The field theoretic Hamiltonian can be written as 

■H^Hi+H2 + H3 + H4 + H5, (29) 

with 

Wi = /dr>t(r")^(r'')V;(r"), (30) 



7^2 = / dr''^\r'')h{r')iP{r''), (31) 



^3 - ^ / / dr?dr^V'^K)^t(ra)^(^a _ r^)^(r^)^(r-) (32) 



^4 = i / / dr\dr'^^Hr\)i>Hr'^)U{r\ - r\)i>{r\)i^{r\) (33) 



rti^ i i dr''dr^^\r'')-^\r^)W{r'' -r^)^{r^)'ij^{r''), (34) 



From the Hamiltonian (29), the equation of motion oi i]:{r°-)i]:{r ) is obtained as 



ifi§M£lMh:lM = [h{r'') + h{r^) + Jdr'''iJj''{r'",t)U{r'' - r'")ijj{r'" ,t) 
+ J dr^'{l}J{r^\ t)U{r^ - r''')'^{r^\ t) 
+ / dr'''Vi^(r"', t)W{r''' - r^)'^{r''\ t) 
+ J dr^'i^^r^' , t)W{r'' - r^')'^{r^\ t)]%jj{r'', t)i;{r\ t), 



(35) 



from which, using ^ dr'^ip'^ (r"" ^t)^{r°' ,t) = ^ dr^ip'^ [r^ ,t)^{r^ ,t) = N, one obtains the equation of motion for the 
joint condensate wavefunction (l){r'^ , r^ , t) , as a generahzation of the Gross-Pitaevskii equation, 

ih^^irlflM = [h(r'') + hir'') + A^//dr°'dr'''|(/.(r'"',r''',t)|2t/(r'^ - r"') 

+NJ J dr'''dr'''\^{r'", r^' , t)\^U{r^ - r^') ,^^. 

+iV//dr'^'dr''V(r'^',r''',t)|2T/F(r°'-r'') ^ ' 

+N f f dr'''dr'''\(t){r'", r''' , t)\^W{r'' - r''')]0(r'^, r^ t). 

This equation may also be obtained directly from the many-body Schrodinger equation, using (p3) and (pq). 

5. QUANTUM COMPUTATION WITH DOUBLE- WELL CONDENSATES 

5.1. Possible robustness 

If a qubit is implemented as a Bose condensate, then there is possible robustness and stability due to the macroscopic 
occupation of a same single particle state. For example, consider an excited state where only one particle is away 
from the ground state, for distinguishable particles there are N possibilities, while there is only one possibility when 
all the particles are identical [g^- For a single condensate, the Bose-condensed state is a product of the same single 
particle state. For two coupled condensates, the Bose-condensed state is now suggested to be a product of the same 
single particle pair state, with symmetrization over all possible ways of pairing. In general, symmetrization always 
needs to be made on the many-body state. This reduces the error probability. Symmetrization has been studied as a 
way of reducing error in quantum computation, which was found to suppress the error probability by 1/A^ [|l3| . Bose 
condensation can be viewed as a natural realization of this prescription. The error reduction due to symmetrization 
may help one to understand why the single particle state emerges out as a macroscopic wavefunction. The robustness 
of condensate wavefunction was demonstrated in the interference experiments [ESl. Nevertheless, the condensate 
wavefunction brings the issue of phase diffusion, which may cause error. 

5.2. NP-complete problems and Nonlinear Quantum Computation 

The class of NP-complete problems is a foundation of the computational complexity theory. It includes thousands 
of practically interesting problems, such as travelling salesman, satisfiability, etc. NP stands for "non-deterministic 
polynomial time" . NP-complete problems are those for which a potential solution can be verified in polynomial time, 
yet finding a solution appears to require exponential time in the worst case. The completeness means that if an 



efficient, i.e. polynomial-time, algorithm could be found for solving one of these problems, one would immediately 
have an efficient algorithm for all NP-complete problems. A fundamental conjecture in classical computation is that 
no such an efficient algorithm exists. Abrams and Lloyd found that with nonlinearity, a quantum computer can 
solve NP-complete problems by efficiently determining if there exists an x for which f{x) = 1, and can solve #P 
problems by efficiently determining the number of solutions [15[. Their algorithm is based on one or two one-bit 
nonlinear gates, together with linear gates. However, it is an experimental fact that fundamental quantum mechanics 
is linear to the available accuracy, while nonlinear fundamental quantum theory p9| usually violates the second law 
of thermodynamics [BOJ and the theory of relativity 13^] . 

Now, because the condensate wavefunction is nonlinear when there is interaction between particles, while it has been 
shown that the condensate wavefunctions of different Bose condensates can be entangled, we propose that condensate 
qubit may be used to realize the nonlinear quantum computing, and thus deal with NP-complete and #P problems. 
It might be also possible to use it to go around some constraints in quantum information processing originated in 
linearity of quantum mechanics. For a Bose condensate of trapped atoms, the atom-atom interaction, and thus the 
nonlinearity of condensate wavefunction, can be tuned. Therefore in principle one may construct linear gates by 
turning off the nonlinearity, and nonlinear gates by turning on the nonlinearity. The nonlinearity of Gross-Pitaevskii 
equation is just of Weinberg- type ||2^, which is used in the algorithms in ||l5| . 

As explained in Sec. 2, in general the condensate wavefunction is not the total pure state of a closed system, hence 
of course its nonlinearity has nothing to do with, and does not share the problems of, nonlinear quantum mechanics. 
In the absence of the interaction, the condensate wavefunction reduces to the pure state of a single particle, and, 
consistently, the nonlinearity disappears. 

5.3. One-bit operation 

We illustrate how to implement the one-bit gates, in terms of a Bose condensate of trapped atoms in a symmetric 
double-well trapping potential V{r) (Fig. 1). This is essentially the Josephson-like effect investigated previously p4[ . 
We may represent bits |0) and |1) as the condensate wavefunctions at the two wells, respectively. Thus 

Mr)\r)dr, (37) 

where n = 0, 1, (f)n{i') — u{r — r„) is the condensate wavefunction corresponding to the local potential V{r — r„), 
which may be parabolic, at the vicinity of the bottom r„. In accordance with Eq. (p^), a qubit \q{t)) is in general a 
superposition of |0) and |1), 

k(t))=co(i)|0)+ci(i)|l)=(^^^|^j), (38) 

where = denotes the matrix representation, with 



io)hJ)ji)-(; 



Gross-Pitaevskii equation leads to 



iH^\q{t))=H\q{t)), (39) 



with 



A.E/+.(;j)+.«(if 1^°,), ,40, 

where 

n = J dr^oirWir) - V{r - ro)]0i(r) (41) 

represents the Josephson-like tunnelling effect, 

7 



K^gJdr\Mr)\\ (42) 

E is the energy of the basis |0) and |1). Therefore, when the s-wave interaction is turned off, k — 0, we may have an 
arbitrary one-bit linear transformation, depending on the time span r: 

\qiT)) ^ e^-^iEI + na,)T]\qiO)). (43) 

Thus one may construct one-bit hnear gates. When the s-wave interaction is turned on, k ^ 0, there is a twisting 
rotation in the state space spanned by |0) and |1). By choosing appropriate time span, this may be used to construct 
one-bit nonlinear gates. 

We mention that when the wavefunction 4>n{i') is real, (t)o{r) ± (t)i{r) are orthogonal, and might be used as qubit. 

5.4. Two-bit Operations 

The nonlinear one-bit gates of Bose-Einstein condensates may be integrated with the linear gates of other qubit 
carriers, so that the algorithms in |15| can be implemented, since only one-bit nonlinear gates are needed there. A 
network of condensates is also possible. In the following, we investigate the evolution of two coupled condensates 
based on a direct interaction, as described in Section 5. We consider trapped atoms in a double- well potential. A long 
range interaction, denoted as W{r — r'), such as dipole-dipole interaction, is a possible basis of the inter-condensate 
interaction. It is considerable either for the magnetic moments of trapped atoms with high magnetic moments [p2| , 
or for the electric dipoles induced by strong dc fields |Q . It would be ideal if the kind of interaction between atoms 
in different condensates, for the purpose of coupling, is absent or somehow canceled out between atoms in a same 
condensate. However, we shall discuss the general case, since it is interesting no matter whether it is used for quantum 
computation. The discussion is formal, without detailed consideration of the suitable physical conditions. 

In the presence of a long range interaction W{r — r') in addition to the s-wave interaction, one should substitute 
U{r — r') — g6{r — r') + W{r — r') in the equation of motion. One-bit Hamiltonian (Eo) is then modified to 

vi oy \^ «;|cip + |ciiVi + |coiV2y ^ ' 

where /^i — jW{r — r')|0„(r)p|(/)„(r')p(ir(ir' is due to W within a same well, while /i2 = jW{r — 
r')\4>n{r)\'^\(j>n{r')\'^drdr' , with n = 1 — n , is due to W within different wells, /zi '^ /X2, since W between atoms 
within a same well is much larger than that between atoms in different wells. 

A two-bit gate may be constructed by putting together two double-wells, each of which confines a condensate. They 
are close to each other in a face-to-face way (Fig. 1), i.e., |0)a is close to |0)6, and |1)q is close to |l)b. Therefore for 
atoms in different condensates, W between, say, an atom in |0)a and an atom in |0)b is much larger than that between 
an atom in |0)a and an atom in |l)f,. This is the origin of the conditional dynamics. 

Substituting Eq. (^^ to the equation of motion of the total condensate wavefunction, Eq. (p6[), with U{r^ — r' ) = 
g^S{r^ — r^ ) + W{r^ — r^ ), {i — a,b), we obtain 

'^^m^ = ^ C„a_„5(/ + II + III + IV + V + VI + VII + VIII), (45) 

with 

/ = / / dr-drV*. (r"^)C. {r'YKr^^ni (r'^)0„5 (r^) ^^^^ 

// = / / dr-dr'rn^ (r")</-*, (r-'')A(r^)^„j (r'^)0„.^ (r^) ^^^^ 

III = g'^N J J J J dr''dr^dr<"dr^'(l)la {r'')(j)l^ {r^)(l)ni (r°)(5f>„5 (r'')|0(r"', r^')\'^5{r'' - r"') 

''5nj.n"<5„5,„.5'^iV|C„.,„.|2 /!</)„. (r-)|4dr-/|</.^.(r'')|4dr-'' (48) 



= 5r,^,n-5^b i>g<'N\C^ 



i^k'^k'' 



IV =NJ J J J dr''dr^dr'''dr^'<j)l^ (^")C'' (^'')</'n? (^")</'n'' (^'')I0(^"', r^')?W{r'' - r"') 



(49) 






2j(rb_rfc') 



(50) 



VI ^N J J J J dr''dr''dr'''dr'>'(j)*„^ (^'')C'. (^'')'/'"? (^")'/'nJ (r'')l'/'(^'", r*"') |2vi^(r'' 



r"') 



(51) 



V^// =N J J J J dr^dr^dr-^'dr^'cljla (r'^)./),*, (r'')(/)„a (r'')<?!)„b (r'')|0(r'^', r''') p VF(r'^' - r'') 



dn" n"d„o , 



A^Eno'd^ 



ICn 



^)- 



n'-" .n^ 1 



(52) 



F/// =N J J J J dr-^dr^dr^'dr^'cljla (r'^m, (r'')0„j (r°)(/'„5 (r*") |<?!)(r'^', r^')\^W{r''' - 
== 5„.,„.5„5,„.A^E„-,„^' |C„»^„,f /l0„.Kr''Vl</>n"(V'^)pM^(r^' - r'^)dr'^cir''' 

= 5„a „a(5„6_„6iV^^(,,(|Co,„t'| +|Cl,„j'| )l'„6'_„a, 



(53) 



where E'', fi*, and k^ 



a, 6) have the same meanings, respectively, as those quantities without the superscripts, 



defined above for a single condensate, n* 



1 



For simplicity, we may set i?" = E , f7" 



r!" 



K" = K". fi^i n'' = / I0n' (■'■') Pl0n''('''* )pW^('''* " ''"* )dr^dr^ is due to the interaction W within a same condensate 
qubit. j/„i „i' = / |0„i(r*)pT4^(r^ — r*)|(/)„j'(r-' )p(ir*dr-' is due to the interaction W between different condensate 
qubits. With symmetry, /i„i „i' = fi\ for n* = n' while fJ.ni,n'' — M2 ^^^^ '^* 7^ ^* • Mi ^ M2- Similarly, i^ni „3' = vi for 
71* = n^ while I'n'.nj' = ^2 for ^' 7^ ^"' • I'l ^ 1^2- Then Eq. (Esi) can be written as a matrix equation. 



where 



ih- 



dt 



fCoo\ 
\CnJ 



( Goo + A^Foo 

n'' G, 




01 + iVi^oi 

Gio + iVi^io 



\ 




n" 

Gn + NFn J 



fCoo\ 

Gqi 

Gio 
\CnJ 



(54) 



G„.„. =E'' + E'' + iff + g'')NK''K'\C^.^,\ 



(55) 



and 



Foo = (/ii + ^i)(2|Goo|' + IGoip + IGioh + (M2 + ^2)(|Goip + |GioP + 2|Giin, 
i^oi = (a*i + ;^2)(|Goo|' + 2\Coif + |Gii|2) + {^l2 + ^^i)(|Goo|' + 2\Ciof + 2\Cnf), 
Foi = (a*i + i^2)(|GooP + 2|Giop + \Cn\^) + (^2 + ^^i)(|GooP + 2|Goip + 2|Giin, 

i^ii = (/ii + ^^i)(|Goip + IGiop + 2|Giin + im + ^2)(2|Goop + |Goip + IGioP). 

M? 



(56) 
(57) 
(58) 
(59) 



We have set /j.^ — ^ — ^i, fi2 = f^2 — f^^- 

In principle, (p4) is a basis for two-bit operations. For the sake of quantum computation, there are some questions 
worthy of investigations, for example, whether (p4) can be used to realize universal two-bit gates, whether there is 
universality for nonlinear gates, how to construct algorithms for NP-complete and #P problems based directly on (|39| ) 
and (^4|), how to realize linear and simpler two-bit operations for the Bose-Einstein condensates, whether swapping 
operation |3J] can be constructed, etc. 



6. SPINOR CONDENSATES 
6.1. Spinor condensate wavefunction 

Up to now, the internal state is irrelevant. In this section, in parallel to the above discussions on a condensate 
in a double-well potential, we discuss the Josephson-coupled internal states of a condensate. In this case, the field 
operator and the condensate wavefunction arc spinors. We use a two-component condensate in an atom trap ||35| as 
the prototype for discussions. 

Suppose that the two internal states are |0) and |1). Then together with the motional degree of freedom, the single 
particle basis state can be written as (/)Q_„(r)|n), (n — 0, 1). The field operator can be defined as 

ip{r,t) = ^ ^ 0Q,n(r')|n)aQ,„(t) 

a n— 0,1 

= Mr,t)\0)+Mr,t)\l), 



y . X ; . (60) 

where 'ipn{r,t) — ^^ 4'a,n{'>')o,a{t). Imposing SGSB on (|6(]| ) leads to the spinor condensate wavefunction 

$(r,t) = a>o(r,i)|0) +<i>i(r,i)|l) = (^ |^|^^Jj ) , (61) 

with <I>(r,i) = (V'(r,t)), <f>„(r,i) = (V'„(r, i)). (|l|) can be written as 

0(r,i) = 0o(r,t)|O)+0i(r,i)|l), (62) 



with 4>{r,t) = ^{r,t)/\/N, 4}nir,t) = <I>„(r, t)/v A^. We have used notations similar to those for the double-well 
condensate. But note the differences in meaning. 

The justification in terms of ODLRO can also be made. The one-particle reduced density matrix is 

(r'|pi|r) = |0)(0|(r',0|pi|r,0) + |0)(l|(r',0|pi|r,l) + |l)(0|(r',l|pi|r,0) + |l)(l|(r',l|pi|r,l), (63) 

where {r',n'\pi\r,n) = (-0„/(r', i)V'J,(r, t)). The existence of ODLRO implies {r',n'\pi\r,n) -^ $„/(r', i)$,* (r, t), and 
thus 

(r'Ipilr) ^ ($o(r', t)|0) -t- <i>i(r', t)\l)){Mr, t)|0) + ^i{r, i)|l))*, (64) 

which indicates the existence of the spinor condensate wavefunction ^(r), as in (pi 



6.2. Four-component spinor condensate wavefunction of twro coupled condensates 

For two coupled two-component Bose systems a and 6, 

^(r^ t)4>{r\ t)^Y.Y.Y.Yl '^""^«" {rn^o.^n^{r^)K) |n'')a„o,„„ (i)a„.,„. {t) 

n" „6 Q," „b 



n" „6 Q," „b 

n") 

^^o(r•^t)^o(r•^t)\ 
^(r^t)V'l(r^t) 

6 



V'l(r•^t)1^o(r•^i) 

V^l(r^i)V'l(r^i)y 



(65) 



Imposing SGSB on (p3q), one obtains 
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$(r^r^t)- Y, E *n^n^(r•^r^i)|n'^) 

n"=0,l„6=o,l 



/$oo(r•^r^t)\ 
$ol(r•^r•^i) 
$lo(r^r^i) 

V$n(r•^r•^t)/ 



(66) 



where ^{r°',r^,t) = {ip{r°-,r^,t)), <&„» „6(r'^,r^, t) = (-(/^na „6(r",r^, i)). Therefore for two coupled two-component 
condensates, the total condensate wavefunction is a four-component spinor. 

The reasoning can also be cast in the form of ODLRO. The one-particle-pair reduced density matrix is 






71°-' 71° ,11°- ,71° 



where 



{r^',n'^\r\n'> |pl|r^ n^r^ n") = {^P^.,ir^' ,t)4,„,,{r' ,tW,{r',t)4>Ur'',t)). 



With ODLRO, 



?,t /„b ,\JA 



(V'„.'K',OV'„-(r'' ^mlir", t)^Ur\t)) ^ <I>„„,^„.,K', r\ t)$:.^„.(r^r^ t) 



Consequently, (q^ approaches 

n°'=0,l r!'''=0,l r!"=04 nt^O,! 

implying the existence of the joint condensate wavefunction as a four component spinor, as given in 
We may write ( pq) as 

</)(r^r^t)= E E '^«^„^(r•^r•^i)|n'^)|r^''), 

n»=0, 1,16=0,1 



(67) 



(68) 



(69) 



(70) 



with <^(r», r^ i) = <I>(r-, r^ t)/N, 0„„,„. (r°, r^ t) = $„.^„. (r», r^ i)/iV. 

The many-body states are similar to the cases in Sec. 3, only with each single particle wavefunction (jja^nir) replaced 
as the combination of motional state and the internal state (j)a,n{r)\n). 



6.3. Hamiltonians and equations of motion 



When there is no coupling between the internal state and the motional state, (j)a,n{i') in ( pO[ ) is independent of n. 
(t)n{i', t) in ( |62| ) can be written as c„(t)4>{r, t), where (/)(r, t) is independent of n, while c„(t) is independent of r. With 
an electromagnetic field, the Hamiltonian for the internal state is 



Hi, 



-CTrr- + -CTz 



(71) 



with o'zln) — (2n — 1)|»t.). uj is the Rabi frequency, while S is the detuning. 

With the coupling between the motional and internal states, the total Hamiltonian for a two-component condensate 
is 



N 



H = El'iC'^O ® i™ + i ® Hra{i)] + E U(^^ - '^i)' 



(72) 



2=1 i<j 

where, with the number of particles in each internal state conserved, the potential U is 



U ■ 



/C/00 \ 

if/01 fC/01 

|f/oi |c/oi 

V UiiJ 



(73) 
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in the basis |00), |01), |10), |11). The field theoretical Hamiltonian is 

H == / rfri/.t(r)[/i(r) (g) Un]i^{r) + / drtP^{r)[i (g> Hin]i^{r 
1 



, , dridr2ip^{ri)ip''{r2)U{ri-r2)il^ir2)ip{ri), (74) 



where the field operator ip{r) is as defined in (|6C|). The equation of motion of ij}(r,t) is 

^^9V^ ^ [(^-iyl + V{r)) ® i.„ + i ® 7?.„ + j i,Hr',t)U{r - r')^(r', t)dr']^ir, t), (75) 

where a 2 x 2 unit matrix is omitted in front of "0^ (r' , t) (similar is the following) . The vertical position of the trapping 
potential depends on the internal state [p5[ , 

V{r) = dtag[Vo{r),Vi{r)], (76) 

where Vn{r) — ujpp'^/2 + uJz{z — z„)^/2. Therefore we have 

,. 5_ / i;o{r,t) \_{Hoo + S/2 C./2 \ / Mr,t) \ 

where 

^„„ = -^ + K(r)+ V j\^n'{r'MUnn'{r-r')dr\ (78) 

2m -^-^ / 

where Uiq = Um. Consequently, through SGSB or ODLRO, one obtains the coupled Gross-Pitaevskii Equations for 
the spinor condensate wavefunction 0„ — {'tpnir)) / ^/N , 

'""dtyMr^t)) { ^/2 Hn-S/2j[Mr,t)J- ^""^ 

Hnn = -^^ + Vn{r)+N ^ j \<j),Ar' ,t)\^Unn' [v ~ r')dr' . (80) 

If the interaction is s-wave interaction, Unn'if — f') = 9nn'S{r — r'). 

For two coupled Bose systems a and 6, the general Hamiltonian is similar to Eq. (28), with each term now changed 
to a matrix in the internal state space, in a way similar to the above case of a single Bose system. In terms of the field 
operators ^p{r°',t) and tp{r^ ^t), formally the field theoretical Hamiltonian can be expressed as (P9|), while the equation 
of motion is in the form of (|3g). However, each potential or interaction energy operator depends in the appropriate 
way on the internal states. In the matrix representation, 

V{r')^dzag[Vo{r'),V,{r')], (81) 

in the basis |0'), |1*), (i = a, fo), while it is assumed that 

W{r' ^r^') = diag[WiiQ{r' ~ri'),WQi{r^ -^r^'^.Wwir^ -r^'),Wii{r' -r^')] (82) 

in the basis of |0*0-'), |0*1-'), ll^O-'), ll'l-'), (j,j = a,b, i ^ j). The long-range interaction within a same condensate 
can be neglected, compared with U. Then it is easy to write Eqs. ( pO| ) to ( p5[ ) in the matrix form. The following is 
the equation of motion for ?/;„a (r", i)?/'„b(r'', t). 



^^ m - y^2^ Vr» +yn-{r )- 2^\/^, +Vn>'{r ) 



+ / dr'^'ijla, (r"', i)[/„o„a' (r" - r-^')^™- {r'" , t) 

+ J dr^'^l (r^', i)C/„.„., (r^ - r^')^„.- (r"', t) (83) 

+ / rfr-" Vt, (r-', t)W^.,,,. (r^ - r-')^n- (r'", t) 

+ / dr^VI.' {r"', t)W„.^., (r- - r^')^„,, (r^', t)]^„. (r", i)^„. (r•^ t). 
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from which one obtains the equation of motion for the condensate wavefunction component (j)^a„b{r°', r ), 

^^ dt i~2;ri Vr" +Vn-{r J - 27;! V^i. +K»(r ) 

+NJJ dr-'dr'''(|0„o,o(r-', r"' , t)\^ + |0„.a(r'^', r^', t)P)t/„.„., (r« - r"') 

+NJJ dr'''dr^'{\ct>ony {r"', r^' M + |0i„.' (r'", r''', t)|2)f/(r^ - r''') (84) 

+iV / / dr-'dr^'(|0o„., (r-', r^', t) ^ + |0,„,, (r-', r"', t)|2)PF(r- - r^')]0„„„. (r", r^ i). 

In principle, the internal state may encode qubit. A simple implementation of one-bit linear operations can be made 
if there is no coupling with the motional state. Then applying an electromagnetic field can realize one-bit operations, 
based on Eq. (|7l|). There is no nonlinearity here, since nonlinearity appears only when the motional degree of freedom 
is involved. 

Because for different bases internal states of a condensate in a trap, vertical positions are different, the inter- 
condensate interaction W,ia„6(r'' — r^) depends on n" and n^ . It can be arranged in such a way (Fig. 2) that 
W^oo('"'^ — '"'') = Wii(r'^ — r**), Woi(r'' — r^) = 1^01(7"° — r**). Then conditional dynamics described above may be 
realized. However, it is difficult to use this to realize two-bit operations for the purpose of quantum computation. One 
might be tempted to realize conditional phases coming from the motional degree of freedom, in a way similar to some 
proposals for individual atoms and ions |^,^. However, the nonlinearity makes this scheme hard to be realized. 

7. Additional remarks 

[1] The effect of nonlinearity on the interference [Q needs more investigations, both from the perspective of physics 
and the perspective of quantum computation. In our scheme of entangling condensates, nonlinearity appears in the 
two-bit evolution. One might find an alternative way, which does not introduce two-bit nonlinearity. On the other 
hand, it is interesting to study computational issues in presence of nonlinearity. Even though two-bit nonlinearity is 
a disadvantage for quantum computing, the physics problem itself is still interesting. 



[2] In the previously studied two-species condensate |39|, each condensate has a single condensate wavefunction, so 



the issue of entanglement is out of question. Moreover, if the two species are two different internal states of a same kind 
atom, they can be Josephson-coupled. A possible way of realizing a mixture with entanglement or a four-component 
spinor condensate, as discussed here, is to mix two different kinds of atoms, each with the two internal states coupled 
through an electromagnetic field. The equation of motion of the joint field can still be written as Eq. (p3[), but 
the main interactions are all due to s-wave scattering, hence W^„a„i. = [/„»„(>, C/„i„3'(r* — r^) = g„in3'S{r^ — r^), 
{ij = a,b). 

[3] An extreme case of the coupled many-particle systems considered here is a system of identical composite particles. 
Then the proposed situation of quantum computation becomes Bose condensation of identical quantum computers. 

[4] Superposition of condensate wavefunctions might be useful for quantum computation no matter whether the 
entanglement is realized and used. 

[5] If a swapping operation can be realized, then probably by using an optical lattice trapping many condensates, 
an architect similar to the one in [p7[ may be constructed, where a head qubit mediates operations between distant 
qubits. 

[6] Many studies on Josephson junction between superconductors are based on quantizing a macroscopic Hamiltonian 
with macroscopic variables such as the charge or particle number. This is equivalent to the approach of condensate 
wavefunction I40L Therefore it is also a kind of condensate qubit that is used in the Josephson-junction quantum 



computation |l4l[42|. Nevertheless, in these proposals, it is either the particle number or the phase of the condensate 
wavefunction that encodes the qubit. In our proposal, it is the two branches of the condensate wavefunction that 
encode the bits. Josephson-junction qubits, since they are charged, were proposed to be coupled via an electric circuit. 

8. SUMMARY 

Coherent properties of Bose-Einstein condensation are well described by condensate wavefunction, which, in a mean 
field theory, is the single particle state in which the condensation occurs. The many-body state is the product of this 
same state occupied by all the particles. In Josephson-like effect, the condensate wavefunction is the superposition 
of two bases wavefunctions. Here we go a step further to suggest that for two such condensates interacting with each 
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other, there is a joint condensate wavefunction which is a superposition of the products of the bases wavefunctions of 
the two condensate, hence there can be entanglement. The many-body state is suggested to be a product of copies of 
this joint condensate wavefunction, with symmetrization. 

With superposition and entanglement, condensate wavefunction, or macroscopic wavefunction, may be used to 
implement quantum computation. That many identical particles occupy a same single particle state may lead to 
some intrinsic robustness and stability, although there is an issue of phase diffusion of the condensate wavefunction. 
For better accuracy, one may try to realize linear evolution of the the condensate wavefunction. On the other 
hand, the nonlinearity due to particle-particle interaction may turn out to be a resource of computational power, as 
indicated in a nonlinear quantum algorithm for NP-complete problems flq]. We have illustrated the ideas by using 
Bose condensation of trapped atoms, especially in doublc-wcU potentials. We have also discussed the existence of a 
four-component condensate wavefunction, as a result of the coupling between two two-component condensates, each 
of which has two Josephson-coupled internal states. 

With these justifications, future researches may include detailed calculations under realistic physical conditions and 
the open issues mentioned above. It is also interesting to extend the consideration to other macroscopic quantum 
coherent systems. It seems that our work in the meantime is the first discussion on quantum entanglement between 
"second quantized" many-particle systems or quantum fields, with off-diagonal long-range order or spontaneous sym- 
metry breaking. Through this work, it is also seen that new physics emergent on a new level of complexity may lead 
to new properties of computation. Indeed, "more is different" |{4^ . 
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APPENDIX 



Here we derive Gross-Pitaevskii equation from ODLRO. 
The one-particle reduced density matrix can be written as 



r'i\piit)\ri) = ,^_ ^., / ■■• / c?^2---dr-Ar(r'i,r2,---,rjv|p(t)|ri,r2,---,rjv), (85) 



(7V-1) 
from which we obtain 



^f^-g^{r'l\Pl{t)\rl} = _ ■■■ dr2- ■ ■ drN{r\,r2,- ■ ■ ,rM\{Hp - pH)\ri,r2,- ■ ■ .tm)- (86) 



(iV-1)! 
One can obtain the equation of motion of pi as 

^h^^{r\\p,{t)\r{) = h{r\){r\\p,{t)\r{) - {r',\p,{t)\r,)h{r,) 

+ j dr2[V{r'^-r2){r[,r2\p2\ri,r2)-{r[,r2\p2\ri,r2)V{r^-r2)], (87) 

where (r']^,r2|/52|ri, r2) is the two-particle reduced density matrix, to which the major contribution comes from 
(ri|pi|ri)(r2|/5i|r2). Assuming 

(r;,r2|p2|ri,r2) = (r'i|/5i|ri)(r2|pi|r2), (88) 

and considering V{ri — r2) — g5{ri — r2), and ODLRO as given in equation M), one obtains 

+gN\^{r[,t)\^^{r[,t)r{ri,t) - gN\4>{r,,t)\^rirum{r[,t), (89) 

which leads to the Gross-Pitaevskii equation. 
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Figure captions: 

Fig. 1. Two interacting Bose condensates. Each condensate is trapped in a double- well potential, hence may 
represent a qubit. |0) is represented by the condensate wavefunction at one well, while |1) is represented by the 
condensate wavefunction at the other well. A two-bit operation between qubits a and b may be based on an interaction 
between atoms in different condensates. 

Fig2. Two interacting two-component Bose condensates. Each condensate has two Joscphson-coupled internal 
states. The trapping potentials for |0) and |1) are displaced between each other, hence the interaction between a 
and b, based on the long-range interaction between atoms in different condensates, depends on the internal states of 
a and b. Each internal state may be regarded as a bit. However, due to coupling with the motional states and the 
nonlincarity, it is difficult to use this situation to realize a two-bit phase gate. 
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